Dissecting Complex Traits Using the Drosophila Synthetic Population Resource

For most complex traits we have a poor understanding of the positions, phenotypic effects, and population frequencies of the underlying genetic variants contributing to their variation. Recently, several groups have developed multi-parent advanced intercross mapping panels in different model organisms in an attempt to improve our ability to characterize causative genetic variants. These panels are powerful and are particularly well suited to the dissection of phenotypic variation generated by rare alleles and loci segregating multiple functional alleles. We describe studies using one such panel, the Drosophila Synthetic Population Resource, and the implications for our understanding of the genetic basis of complex traits. In particular, we note that many loci of large effect appear to be multiallelic. If multiallelism is a general rule, analytical approaches designed to identify multiallelic variants should be a priority for both genome wide association studies and multi-parental panels

How Good a Deal Was the Tobacco Settlement?: Assessing Payments to Massachusetts

We estimate the increment in Massachusetts Medicaid program costs attributable to smoking from December 20, 1991, to 1998. We describe how our methods improve upon earlier estimates of analogous costs at the national level. Current costs to the Massachusetts Medicaid program approximate the payments to Massachusetts under the tobacco settlement of November 1998. Whether these payments are viewed as appropriate compensation for Medicaid costs over time depends upon the rate of increase in future health care costs, the rate of decline in smoking, the proportion of smoking that should be attributed to the actions of the tobacco companies and the liklihood that state would have prevailed at trial. The costs to the Medicaid program are dwarfed by the internal costs to smokers themselves.

The Genetic Architecture of Methotrexate Toxicity Is Similar in Drosophila melanogaster and Humans

The severity of the toxic side effects of chemotherapy varies among patients, and much of this variation is likely genetically based. Here, we use the model system Drosophila melanogaster to genetically dissect the toxicity of methotrexate (MTX), a drug used primarily to treat childhood acute lymphoblastic leukemia and rheumatoid arthritis. We use the Drosophila Synthetic Population Resource, a panel of recombinant inbred lines derived from a multiparent advanced intercross, and quantify MTX toxicity as a reduction in female fecundity. We identify three quantitative trait loci (QTL) affecting MTX toxicity; two colocalize with the fly orthologs of human genes believed to mediate MTX toxicity and one is a novel MTX toxicity gene with a human ortholog. A fourth suggestive QTL spans a centromere. Local single-marker association scans of candidate gene exons fail to implicate amino acid variants as the causative single-nucleotide polymorphisms, and we therefore hypothesize the causative variation is regulatory. In addition, the effects at our mapped QTL do not conform to a simple biallelic pattern, suggesting multiple causative factors underlie the QTL mapping results. Consistent with this observation, no single single-nucleotide polymorphism located in or near a candidate gene can explain the QTL mapping signal. Overall, our results validate D. melanogaster as a model for uncovering the genetic basis of chemotoxicity and suggest the genetic basis of MTX toxicity is due to a handful of genes each harboring multiple segregating regulatory factors

Functional Expression of Parasite Drug Targets and Their Human Orthologs in Yeast

BACKGROUND: The exacting nutritional requirements and complicated life cycles of parasites mean that they are not always amenable to high-throughput drug screening using automated procedures. Therefore, we have engineered the yeast Saccharomyces cerevisiae to act as a surrogate for expressing anti-parasitic targets from a range of biomedically important pathogens, to facilitate the rapid identification of new therapeutic agents. METHODOLOGY/PRINCIPAL FINDINGS: Using pyrimethamine/dihydrofolate reductase (DHFR) as a model parasite drug/drug target system, we explore the potential of engineered yeast strains (expressing DHFR enzymes from Plasmodium falciparum, P. vivax, Homo sapiens, Schistosoma mansoni, Leishmania major, Trypanosoma brucei and T. cruzi) to exhibit appropriate differential sensitivity to pyrimethamine. Here, we demonstrate that yeast strains (lacking the major drug efflux pump, Pdr5p) expressing yeast ((Sc)DFR1), human ((Hs)DHFR), Schistosoma ((Sm)DHFR), and Trypanosoma ((Tb)DHFR and (Tc)DHFR) DHFRs are insensitive to pyrimethamine treatment, whereas yeast strains producing Plasmodium ((Pf)DHFR and (Pv)DHFR) DHFRs are hypersensitive. Reassuringly, yeast strains expressing field-verified, drug-resistant mutants of P. falciparum DHFR ((Pf)dhfr (51I,59R,108N)) are completely insensitive to pyrimethamine, further validating our approach to drug screening. We further show the versatility of the approach by replacing yeast essential genes with other potential drug targets, namely phosphoglycerate kinases (PGKs) and N-myristoyl transferases (NMTs). CONCLUSIONS/SIGNIFICANCE: We have generated a number of yeast strains that can be successfully harnessed for the rapid and selective identification of urgently needed anti-parasitic agents

Genetic Dissection of the Drosophila melanogaster Female Head Transcriptome Reveals Widespread Allelic Heterogeneity

Modern genetic mapping is plagued by the “missing heritability” problem, which refers to the discordance between the estimated heritabilities of quantitative traits and the variance accounted for by mapped causative variants. One major potential explanation for the missing heritability is allelic heterogeneity, in which there are multiple causative variants at each causative gene with only a fraction having been identified. The majority of genome-wide association studies (GWAS) implicitly assume that a single SNP can explain all the variance for a causative locus. However, if allelic heterogeneity is prevalent, a substantial amount of genetic variance will remain unexplained. In this paper, we take a haplotype-based mapping approach and quantify the number of alleles segregating at each locus using a large set of 7922 eQTL contributing to regulatory variation in the Drosophila melanogaster female head. Not only does this study provide a comprehensive eQTL map for a major community genetic resource, the Drosophila Synthetic Population Resource, but it also provides a direct test of the allelic heterogeneity hypothesis. We find that 95% of cis-eQTLs and 78% of trans-eQTLs are due to multiple alleles, demonstrating that allelic heterogeneity is widespread in Drosophila eQTL. Allelic heterogeneity likely contributes significantly to the missing heritability problem common in GWAS studies.This work was supported by NIH R01 RR024862/OD010974 to SJM and ADL, an American Recovery and Reinvestment Act (ARRA) administrative supplement to this award, and F32 GM099382 to EGK. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript

An R Framework for the Partitioning of Linkage Disequilibrium between and Within Populations

Patterns of linkage disequilibrium (LD) across the genome result from a myriad of contributing factors including selection and genetic drift. Natural selection can increase LD near individually selected loci, or it can influence LD between epistatically selected groups of loci. Statistics have previously been derived which compare levels of linkage disequilibrium in subpopulations relative to the total population. These statistics may be leveraged to identify loci that may be under selection or epistatic selection. This is a powerful approach, but to date no framework exists to support its use on a genome-wide scale. We present ohtadstats, an R package designed to facilitate the implementation of Ohta’s D statistics in a variety of use cases. Statistics calculated by this package can be used to determine whether a locus is under selection or not, and can provide insight into the nature of the selection that is taking place (hard sweep or epistatic selection). This package is available on the Comprehensive R Archive Network (CRAN).   Funding statement: This research was supported by funding from the USDA Agricultural Research Service. PFP is funded by the University of Missouri Life Sciences Fellowship and a training grant from the National Institute of Health (T32GM008396)

Structure of the complex of an Fab fragment of a neutralizing antibody with foot-and-mouth disease virus: Positioning of a highly mobile antigenic loop

Data from cryo-electron microscopy and X-ray crystallography have been combined to study the interactions of foot-and-mouth disease virus serotype C (FMDV-C) with a strongly neutralizing monoclonal antibody (mAb) SD6. The mAb SD6 binds to the long flexible GH-loop of viral protein 1 (VP1) which also binds to an integrin receptor. The structure of the virus-Fab complex was determined to 30 Å resolution using cryo-electron microscopy and image analysis. The known structure of FMDV-C, and of the SD6 Fab co-crystallized with a synthetic peptide corresponding to the GH-loop of VP1, were fitted to the cryo-electron microscope density map. The SD6 Fab is seen to project almost radially from the viral surface in an orientation which is only compatible with monovalent binding of the mAb. Even taking into account the mAb hinge and elbow flexibility, it is not possible to model bivalent binding without severely distorting the Fabs. The bound GH-loop is essentially in what has previously been termed the 'up' position in the best fit Fab orientation. The SD6 Fab interacts almost exclusively with the GH-loop of VP1, making very few other contacts with the viral capsid. The position and orientation of the SD6 Fab bound to FMDV-C is in accord with previous immunogenic data.Peer Reviewe

Additional Evaluation of the Point-of-Contact Circulating Cathodic Antigen Assay for Schistosoma mansoni Infection.

Studies of the urine-based point-of-contact cathodic circulating antigen test (POC-CCA) in Schistosoma mansoni-endemic settings in Africa indicate it has good sensitivity in detecting infections, but in areas of low prevalence, the POC-CCA can be positive for persons who are egg-negative by Kato-Katz stool assays. We examined the POC-CCA assay for: (a) batch-to-batch stability; (b) intra-reader and inter-reader variability; (c) day-to-day variability compared to Kato-Katz stool assays, and (d) to see if praziquantel (PZQ) treatment converted Kato-Katz-negative/POC-CCA positive individuals to POC-CCA negativity. We found essentially no batch-to-batch variation, negligible intra-reader variability (2%), and substantial agreement for inter-reader reliability. Some day-to-day variation was observed over 5 days of urine collection, but less than the variation in Kato-Katz stool assays over 3 days. To evaluate the effect of treatment on Kato-Katz(-)/POC-CCA(+) children, 149 children in an area of 10-15% prevalence who were Kato-Katz(-) based on 3 stool samples but POC-CCA(+) were enrolled. Seven days after treatment (PZQ 40 mg/kg) samples were again collected and tested. Almost half (47%) POC-CCA positive children turned negative. Those still POC-CCA positive received a second treatment, and 34% of them turned POC-CCA negative upon this second treatment. Most who remained POC-CCA positive shifted each time to a "lesser" POC-CCA "level of positivity." The data suggest that most Kato-Katz-negative/POC-CCA positive individuals harbor low-intensity infections, and each treatment kills all or some of their adult worms. The data also suggest that when evaluated by a more sensitive assay, the effective cure rates for PZQ are significantly less than those inferred from fecal testing. These findings have public health significance for the mapping and monitoring of Schistosoma infections and in planning the transition from schistosomiasis morbidity control to elimination of transmission

Genetic dissection of a model complex trait using the Drosophila Synthetic Population Resource

Genetic dissection of complex, polygenic trait variation is a key goal of medical and evolutionary genetics. Attempts to identify genetic variants underlying complex traits have been plagued by low mapping resolution in traditional linkage studies, and an inability to identify variants that cumulatively explain the bulk of standing genetic variation in genome-wide association studies (GWAS). Thus, much of the heritability remains unexplained for most complex traits. Here we describe a novel, freely available resource for the Drosophila community consisting of two sets of recombinant inbred lines (RILs), each derived from an advanced generation cross between a different set of eight highly inbred, completely resequenced founders. The Drosophila Synthetic Population Resource (DSPR) has been designed to combine the high mapping resolution offered by multiple generations of recombination, with the high statistical power afforded by a linkage-based design. Here, we detail the properties of the mapping panel of >1600 genotyped RILs, and provide an empirical demonstration of the utility of the approach by genetically dissecting alcohol dehydrogenase (ADH) enzyme activity. We confirm that a large fraction of the variation in this classic quantitative trait is due to allelic variation at the Adh locus, and additionally identify several previously unknown modest-effect trans-acting QTL (quantitative trait loci). Using a unique property of multiparental linkage mapping designs, for each QTL we highlight a relatively small set of candidate causative variants for follow-up work. The DSPR represents an important step toward the ultimate goal of a complete understanding of the genetics of complex traits in the Drosophila model system.This work was supported by the following NIH R01 grants: RR024862 to S.J.M. and A.D.L., GM085260 to S.J.M., GM085251 to A.D.L., GM078338 to S.S., and GM074244 to K.W.B